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In this paper we develop a Markov Chain Monte Carlo code to study the dark 
matter properties in frameworks to interpret the recent observations of cosmic ray 
electron/positron excesses. We assume that the dark matter particles couple domi- 
nantly to leptons and consider two cases, annihilating or decaying into lepton pairs, 
respectively. The constraint on the central density profile from H.E.S.S. observation 
of diffuse 7-rays around the Galactic center is also included in the Markov Chain 
Monte Carlo code self-consistently. In the numerical study, we have considered two 
cases of the background: fixed e'^e~ backgrond and the relaxed one. Two data sets 
of electrons/positrons, i.e. PAMELA+ATIC (Data set I) and PAMELA+Fermi- 
LAT+H.E.S.S. (Data set II), are fitted independently, considering the current incon- 
sistence between the observational data. We find that for the Data set I, dark matter 
with ~ 0.70 TeV for annihilation (or 1.4 TeV for decay) and a non-negligible 
branching ratio to e"^e~ channel is favored; while for the Data set II, ~ 2.2 
TeV for annihilation (or 4.5 TeV for decay) and the combination of ^~ and t'^t~ 
final states can best fit the data. We also show that the background of electrons 
and positrons actually will significantly affect the branching ratios. The H.E.S.S. 
observation of 7-rays in Galactic center ridge puts a strong constraint on the central 
density profile of the dark matter halo for the annihilation dark matter scenario. In 
this case the NFW profile which is regarded as the typical predication from the cold 
dark matter scenario, is excluded with a high significance (> 2>a). For the decaying 
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dark matter scenario, the constraint is much weaker. 
PACS numbers: 98.80.Es; 98.80.Cq 

I. INTRODUCTION 

The recent reported results on the abnormal excesses o: 



by PAMELA and the spectra of electrons^ by ATIC 3, PPB-BETS Q, H.E.S.S. Q, Q 



cosmic ray (CR}_positron fraction 



and Fermi-LAT [6l have invoked extensive discussions on the possible existence of dark 

nh 

matter (DM) signals (for a recent review, see e.g., [7|, |8[). Meanwhile, the ratio of antiproton 
to proton measured by PAMELA 9| is well consistent with the astrophysical expectation 
from interactions between CR nuclei and interstellar medium (ISM) [lO|. It indicates that 
if the DM contributes to the electron/positron excesses, it should dominantly annihilate or 
decay into leptons instead of gauge bosons or quarks [ll|, [3] • Although the observational 
data of the electron spectra from ATIC and Fermi-LAT/H.E.S.S. are not fully consistent, it 
has been shown in literature that the DM models with annihilation or decay modes directly 
to leptons can give good description to the observational data, given proper mass of DM 
particle and flavors of the final state particles (e.g., 13-16t). Specifically, the ATIC data 



favor a e^e channel to describe the bump and fast drop for energy ~ 600 GeV 



While for Fermi-LAT and H.E.S.S. data, a softer electron/positron spectrum from the decay 
of /i^/i" or r'^r" final states can better reproduce the smooth behavior [isl. 

However, for the numerical studies of fitting to the data in the literature so far one usually 
takes some specific parameters for a given model, then fits to the data for illustrations of 
how the model works instead of performing a global analysis. Therefore bias exists in the 
conclusions drawn from this kind of studies. A global fit to the observational data thus will 
be very useful to extract the model-independent implication from the data and to explore 
the correlations among different parameters. In this work we employ a Markov Chain Monte 



Carlo (MCMC) technique |l8l-l20| to globally fit the parameters of the DM scenario to the 
experimental data on the electron/positron spectrum observed. The constraint from the 
diffuse 7-ray emission around the Galactic center (GC), known as GC ridge, by H.E.S.S. 



^ When we mention about experimental data of electron spectra, it actually means the total spectra of 
electrons and positrons. 
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2l| is also included in the MCMC program. Based on the global fitting results, we can then 
further investigate the possible implication on the models of both on the particle physical 
side and also the astrophysical side of DM, in a manner of self-consistency. 

Another issue which is not seriously taken into account in previous studies is the influence 
of electron/positron background on the DM models. The global fit method makes it possible 
to include the uncertainties of the background contribution. This point is also discussed in 
this work. 

The paper is organized as follows. In Sec. II we introduce the production mechanism of 
the electrons and positrons generated from DM annihilation or decay, and their propagation 
effect in the Milky Way (MW). In Sec. Ill we describe the accompanied 7-rays which are 
served as a cross check of the self-consistency of the model configuration. The MCMC fitting 
procedure and results are given in Sec. IV. Finally, Sec. V is the summary. 

II. PRODUCTION AND PROPAGATION OF THE COSMIC e± 

Due to the constraints from p/p by PAMELA, the DM is thought to be dominately 
coupled with leptons [ll], 12|. Since in this work we mainly focus on the methodology of 



using MCMC to globally fit the observational data, we will limit our discussion in the cases 
that DM annihilates or decays to lepton pairs directly^. The production rate of electrons 
(positrons) can be written as 

for DM annihilation, or 

/ ^^ 1 dA^ 

qe{r,E) = — 

m^T ah 

for DM decay, where is the mass of DM particle, (erf) or r are the annihilation cross 
section or decay age of DM respectively, is the electron yield spectrum for one anni- 
hilation or decay with i = e, fi, r, and p(r) is DM spatial density in the MW halo. The 
density profile of the MW halo is taken as the form 

^^"^^ " (r/r,)7(l+r/r,)3-7' 



P(r), (2) 



^ A more general discussion including the quarks or gauge bosons in the final states will be a straightforward 
extension of the current study and will be present in a future publication. 
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where 7 represents the central cusp slope of the density profile, and ps are scale radius 
and density respectively. For the MW DM halo, we adopt the total mass to be Mmw ~ 10^^ 

Mo 



n 



221 ] and the concentration parameter to be Cmw ~ 13.5 [2^^]. Then we have 



'"mw/cmw(2 — 7) where rMw ~ 260 kpc is the virial radius of MW halo. Then ps can be 
derived by requiring Mmw = / pdV. The local density in this process is checked to be 
within 0.27 to 0.25 GeV cm^'^ for 7 varying from to 1.5. 

After the production from DM annihilation or decay, the electrons and positrons will 
propagate diffusively in the MW due to the scattering with the random magnetic field. 
Besides the diffusion, the dominant process of electron propagation is the energy loss from 
synchrotron radiation in Galactic magnetic field and the inverse Compton (IC) scatterings in 
the interstellar radiation field (ISRF). There may also be global convection driven by stellar 



2^, these 



wind and reacceleration by random interstellar shock, however, as shown in Ref. 
effects can be safely neglected for electron energy > 10 GeV. The propagation equation is 



-DV^N{r,E) 



_d_ 
dE 



dt 



■N{r,E) 



q{r,E), 



(4) 



where D{E) = I3DqIZ^ {IZ = pc/Ze is the rigidity) is the diffusion coefficient, dE/dt = 
— c^/te with e = E/1 GeV and te ~ 10^^ s, is the energy loss rate for typical values of the 
Galactic magnetic field and ISRF, g(r, E) is the source function of as given in Eqs. ([1]) 
and ([2]), and N{r,E) is the propagated e"^ spectrum. 

The propagator for a point source located at (r, z) from the solar location with monochro- 
matic injection energy Es can be written as 25|, |26| 



ge{r,z,E^Es] 
in which we define a pseudo time f as 

f = te- 



— X g0(r,z,r). 



(5) 



1-6 



(6) 



^©(r, z, f) is the Green's function for the re-arranged diffusion equation with respect to the 
pseudo time f 



exp 



(7) 



The effect of boundaries along z = ±L appears in only. Following Ref. |25j we use two 
distinct regimes to approach Q^^: 
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for ( = L'^/ADqt ^ 1 (the extension of electron sphere A = \^4Dqt is small) 



o-i.r>^i:i-rr^..,[-^^). (s) 



where Zn = 2Ln + {—l)"'z; 
otherwise 



^ oo 



where 



n=l 



Mz) = sm[kn{L - \z\)]; kn = {n-l/2)Ti/L, (10) 
0;(^) = sin[A;;(L-^)]; k'^ = mi/L. (11) 



For any source function g(r, z, ^; E^) the local observed flux of positrons can be written as 



^0 = — X 2 / / rdr / dEse;0(r, z,E ^ Es) / d^g(r, ^, 6; Eg). (12) 
7o -/o Je Jo 



An 



For the propagation parameters, we use the medium (referred as "MED") set o: 
eters which is derived through fitting the observational B/C data given in Ref. 
Do = 0.0112 kpc^ Myr-\ 5 = 0.70 and the height of the diffusive halo L = 4 kpc. 



param- 



27|, i.e., 



III. GAMMA RAYS 



For the purely leptonic annihilation or decay of DM particles, 7-rays can be generally 
produced in two ways. One is the final state radiation (FSR) which is emitted directly from 
the external legs when DM particles annihilates or decays to charged leptons'^. The other 
is the IC scattering photons from ISRF when the electrons and positrons propagate in the 
MW. Compared with the IC photons, the FSR has several advantages. First the spectrum 
of FSR is unique and may give smoking gun diagnostic for DM signal. In addition, the 
FSR does not depend on the astrophysical environment such as the distribution of ISRF. 
Therefore in this work we only consider the FSR. It will simplify the calculation, meanwhile 
the results obtained on the constraints on the model parameters is also conservative. 



^ Note that for the tau channel, the decay of r leptons will produce a large number of neutral pions which 
can then decay into photons. We also include this contribution in the FSR. 
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The photon yield spectrum for e"^ or /i^ channel for S> rrie, can be written as 



28 



29| 



diV 
dx 



log — 1 

IT X \ mf 



X 



(13) 



where a ~ 1/137 is the fine structure constant, i = e, fi. For DM annihilation we have 



s = 4:m'^ and x = Ej/m^; while for DM decay s = m"^ and x = 2E^/m^ 30[. For r 



channel, we adopt the total parameterization including the direct FSR component as shown 
in Eq.( [T3|) and the decay products from the chain r — )■ 7r° — 7 [31 1 



dx 



X 



-1.31 



(6.94x - 4.93x2 - 0.51x^) e 



4.53a; 



(14) 



with the same definition of above. 

The 7-ray flux along a specific direction can be written as 



0(E^,^) = C X W{E^) X J(^) 



47r 







X 



X 

1 dN 



X 



^ Ilos P"^^^^^^' for annihilation 
;i/LO5P(0d/, for decay 



(15) 



where the integral is taken along the line-of-sight (LOS), W{E) and J{ip) represent the 
particle physics factor and the astrophysical factor respectively, Rq = 8.5 kpc is the solar 
system location from GC, and is the local DM density. For the emission from a diffuse 
region with solid angle Af2, we define the average astrophysical factor as 

1 



J. 



An 



j{i/j)dn, 



(16) 



which is fully determined by the parameter 7 in Eq.(l3l). 

In this work we use the diffuse 7-ray observation in the sky region |/| < 0.8° and |6| < 0.3^^ 



around the GC by H.E.S.S. 



21| to constrain the central profile of DM density distribution. 



It should be noted that the H.E.S.S. observation of the GC ridge 7-ray flux is a background 
subtracted one. In Ref. j2l| the emission from 0.8° < \b\ < 1.5°, |/| < 0.8° is taken as the 
background. Therefore the H.E.S.S. reported result is not the total emission of the selected 
sky region, but a lowered one. To compare the calculated result with the data, we calculate 
the 7-ray flux from DM in the H.E.S.S. signal region (|6| < 0.3°, |/| < 0.8°) with subtracting 
the one in the H.E.S.S. background region (0.8° < |6| < 1.5°, |/| < 0.8°). Thus the Jao 
factor adopted for the 7-ray calculation actually means J^fj — in the following. 
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Similar to Ref. [32[, we employ a power- law component = a^E^y ^ to represent the 
astrophysical background contribution to the 7-rays. The FSR from DM is then added to 
the background to fit the observational data. The fit to the 7-ray data is also included in 
the MCMC code so that we can derive a consistent constraint on the DM density profile. 

IV. MCMC METHOD AND RESULTS 
A. Method 

In our study, we perform a global analysis employing the MCMC technique to determine 
the parameters related to DM models. The MCMC sampler is implemented by using the 
Metropolis-Hastings algorithm. It is an efficient procedure for generating samples which 
lies in the fact that its equilibrium distribution corresponds to the likelihood function in 
parameter space. Comparing with a traditional approach where the likelihood function is 
evaluated on a grid of points in parameter space, the MCMC method has the great potential 
in expanding the dimension of the parameter series, since the computational requirements 
of MCMC procedures are insensitive to the dimensionality of the parameter space. 

Fixing the e+e~ backgrond, we have the following parameter space: 

P = (m^, {av) or r, B^, B^, B^, 7, a^, 6-y), (17) 

where is the DM particle mass, {av) is the annihilation cross section, for decaying DM r 
denotes the lifetime, i?e, B^ and B^ are the annihilation (or decay) branching ratios of DM 
particle into e^e~ , fi^fi~ and r^r~ pairs. For the scenario considered in this work, we have 
the constraint condition Be + Bfj_ + B^ = 1. The parameter 7 parametrizes the DM density 
profile and a^, are parameters characterizing astrophysical 7-ray background. So in our 
global fitting procedure 7 free parameters are in general involved. 

With variables in (fTTIl we can calculate the theoretical expectations of electrons and 
positrons produced from DM. However, to compare with the observational data, we still 
need to give the astrophysical background contribution to the electrons and positrons. Since 
the background calculation, especially for positrons, is quite complicated and consumes a 
lot of time, we employ two approaches of the background. Firstly we fix the background 



calculated by GALPROP package 33| based on a conventional diffusion + convection (DC) 
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model 



12 



34j . It is shown that the conventional GALPROP model can give fairly good 



description to the CR data before PAMELA/ATIC/Fermi-LAT/H.E.S.S. (l2j, as well as the 
all-sky diffuse 7-ray data from EGRET (except the "GeV excess", 35|]) and Fermi-LAT at 
intermediate latitudes SgJ. The GALPROP model parameters are the same as Ref. 12 1 
except for the high energy electron injection spectrum we adopt a slightly softer one 2.58. 
For the second approach we include a more complete treatment on this problem. A power 
law function q^,- = ae-E^^"' with two free parameters a^-, b^- is involved to describe the 
injection source of primary electrons. For background positrons which are thought to be 
produced through CRs interact with ISM, we use the locally measured proton and Helium 
spectra as parameterized in Ref. [l^ and an average ISM density ~ 1 cm~^ to calculate 
the positron production rate. For the secondaries production from p — p inelastic collision 
we use the parameterization given in Ref. js^]- The interaction is restricted in a thin disk 
with half height ~ 0.1 kpc. The spatial distribution of positron source is neglected, which 
is demonstrated not to affect the final positron spectrum 2J]. In addition, we employ a 
free factor Ce+ lies between 0.5 and 2 to describe the possible uncertainties about the ISM 
density and the disk height. That is to say, totaly we will have 10 free parameters in this 
approach. 



B. Data sets 



Since the electron spectrum measurements between ATIC and Fermi-LAT/H.E.S.S. still 
have discrepancy, for an unbiased study at the current stage we will try to fit two differ- 
ent data sets respectively. One is to combine the PAMELA positron fraction data, ATIC 
electron data^ and the H.E.S.S. 7-ray data (Data set I), and the other is the combination 
of PAMELA, Fermi-LAT, H.E.S.S. electron and H.E.S.S. 7-ray data (Data set II). Note 
that for PAMELA we use only the data with energies higher than ~ 5 GeV in the MCMC 
fitting code. The low energy data are thought to be more easily to be affected by the solar 
modulation. We get the by comparing the spectrum between the theoretical expectations 

4 The ATIC data include ATIC1+ATIC2+ATIC4, which are taken from the talk by J. Isbert in TANGO 
in PARIS workshop, see http://irfu.cea.fr/Mcetings/tangoinparis/main.htm, From the numerical calcu- 
lation, we find that the fitting results between ATIC1+ATIC2+ATIC4 and the published ATIC1+ATIC2 
data are very similar, but the former one converges much faster and better. 
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and the corresponding observational values. For the PAMELA, the observation quantity is 
the positron fraction /{e^ + e~), while for ATIC, Fermi-LAT and H.E.S.S. the undistin- 
guished fluxes of positrons and electrons are concerned. We have taken the total likelihood 
to be the products of the separate likelihoods C of all the data. In other words, defining 
X- = -21og£i, we get 

Xtotal ~ ^ ] Xi 1 (18) 
i 

where i labels different observational data. 



C. Numerical Results 

1. Fixed e^e~ background 

Firstly we consider the fixed background approach. The one dimensional probability 
distributions of the fitting parameters for annihilation DM scenario are shown in Fig{Tl The 
left eight panels are for Data set I, and the right eight panels are for Data set II respectively. 
To get the ID probability of each parameter we have marginalized over the other parameters. 
The details about the parameters are compiled in Table HI 

To reproduce the peak of ATIC observation of Data set I, the mass of DM is found to 
be ~ 0.70 TeV with annihilation channel almost purely to e^e^. The 2cr upper limits for 
and r+r" channels are < 0.269 and Bj. < 0.226 respectively. While for the fit 
to Data set II, a larger mass ^ 2.2 TeV with larger B^ and Bj. which can give softer 
electron/positron spectra is favored. The e~^e~ channel which can give strongly peaked 
electron/positron spectra is remarkably suppressed. The 2a upper limit for Be is ~ 0.03. 

An important issue is the constraint from GC diffuse 7-rays observed by H.E.S.S.. A 
strong constraint on the central DM density profile is shown. At 2cr level, we find 7 < 0.539 
for Data set I and < 0.515 for Data set II. It is significantly flatter than the canonical NEW 



profile with 7 



4l|. This indicates that DM particle might not be so cold, but behaves 



like warm DM instead 



32 



|42j . We can also note that the probability distributions of -B^ and 
Br are relatively broad for the Data set II fit. The reason for this is due to the degeneracy 
between and r+r" channels. We will discuss this issue later. 

The results for decaying DM model are shown in Fig. [2]and Table HTl respectively. We can 
see that the results are very similar to the annihilation case except for the constraint on the 
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FIG. 1: Probability distributions of the eight parameters in annihilation DM scenario used to fit 
Data set I (left) and II (right) respectively, for fixed e^e~ approach. To get the ID probability of 
each parameter we have marginalized over the other parameters. 




FIG. 2: The same as Fig{T]but for decaying DM scenario. 



density profile parameter 7. It is just as expected that since the electron/positron source 
function is proportional to p for DM decay instead of the dependence of DM annihilation, 
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TABLE I: Mean 1 a errors or 95% limits for the parameters in annihilation DM scenario for the 
fixed e"'"e~ background approach. 



Parameters 


Data Set I 


Data Set II 


my (TeV) 

A. V / 


0.711+Hi 

— U.U4z 


2.211+°{}? 

— U. ill 


(av) (10"23 cm^ s^M 


0.969+°-j^^ 


9.105+j-j^I 


Be 


0.795 ± 0.090 


< 0.031 




< 0.269 


U.Dyi_Q Q5Q 


Br 


< 0.226 




7 


< 0.539 


< 0.527 


(10-5 GeV-i cm'^ s^^ sy~^) 




< 19.608 






9 9Qft+0.048 
Z.ZOU_Q Q35 



TABLE II: Mean 1 a errors or 95% limits for the parameters in decaying DM scenario for the fixed 
e^e~ background approach. 



Parameters 


Data Set I 


Data Set II 


(TeV) 


i-4i8„o.088 


4 47c;+0.231 


r (10^6 s) 


9 Qc:q+0.252 
^••Jjy_0.246 


n 704+0.066 


Be 


7+0.049 


< 0.025 


B, 


< 0.387 


n ^fi^+O-032 


Br 


< 0.333 


n 421+0-046 
U.4Zi_Q 033 


7 


< 1.114 


< 1.130 


(10^5 GeV-i cm-^ g-i st~^) 


9.936l^iii° 


< 18.498 




97t:+0.030 
Z.Z /O_o.o35 


91 c+0.066 
Z.ZiO_o.o33 



the constraint on the central cusp slope is expected to be much weaker. 

In Fig. [3]we show the results of the total electron + positron spectra and positron fraction 
for the best-fitting parameters, for annihilation DM scenario. The results show a very good 
agreement with the observational data. The results for decaying DM scenario, which are not 
shown here, are almost undistinguishable from the annihilation case. Because the final states 
from DM annihilation and decay are completely the same, the only difference comes from 
the spatial distribution of DM induced electrons/positrons. While the observed high energy 
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electrons/positrons should mainly come from the local regions near the Earth where there 
is no big difference between the distribution for DM annihilation and the p distribution 
for DM decay, therefore the propagated electron/positron spectra are very similar for both. 
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FIG. 3: Left: the e+ + e spectrum inc 
with the observational data from ATIC 



uding the contribution from DM annihilation compared 
, HESS P,3 and Fermi-LAT \ Right: the e+/(e-+e+) 



ratio including the contribution from DM annihilation as a function of energy compared with the 



data from AMS 



HEAT 



3g 



401] and PAMELA 



Finally we investigate the correlation among the various annihilation or decay channels. 
We find that there is no evident correlation between parameter B^. and other two branching 
ratios. This can be understood that the spectrum of electrons/positrons from e^e~ channel 
has a very hard and spiky signature, which can be easily distinguished from the other two 
channels. However, the cases will be different for fi~^p~ and t~^t~ channels. For the p^ p~ 
or T^T~ channels, the electrons are generated through the decay of muons or tauons, which 
both can give a softer electron spectrum. Thus the results from fi~^p~ channel and t~^t~ 
channel should have some degeneracies. In Fig. |4] we plot the 2D correlation between 5^ 
and for the fit to Data set II, for DM annihilation and decay respectively. It indeed 
shows a strong anti-correlation between 5^ and Br- That means to describe the Data set 
II the muon channel is equivalent to the tauon channel to some extent. It will be not easy 
to distinguish whether the DM particle couples mainly with muon flavor or tauon flavor 
through the electron/positron data. This degeneracy will also lead to an uncertainty on 
the study of 7-ray prediction or constraint. For tauon channel, more 7-ray photons can be 
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generated through tauon decay, while for muon channel the FSR radiated 7-ray photons will 
be much fewer. On the other hand, we expect future accurate 7-ray experiment to break 
this degeneracy. 



0.9 
0.8 
0.7 
0.6 

ca" 0.5 
0.4 
0.3 
0.2 
0.1 


( 

FIG. 4: Two dimensional plot for the correlation between parameters and for the fit to Data 
set II, for DM annihilation (solid) and decay (dashed) respectively. The inner and outer contours 
show the results at la and 2a confidence levels. 

2. Varying e~^e~ background 

In this section, we present the results by setting the e~^e~ background free. Three addi- 
tional parameters a^,-, fee- and Ce+ are added in the MCMC approach. The ID parameter 
distributions are shown in Figs. El M and Tables IIIIl IIVI respectively, for the four cases as 
considered in the fixed e~^e~ background section. 

We find that there are remarkable effects on the branching ratios if changing the back- 
ground. This is reasonable since a different background will lead to a different blank for 
DM to fill. For example if the background electron spectrum is harder, then the required 
contribution from DM is expected to be softer. Actually if we relax the background electron 
spectrum, a harder spectrum compared with the previous fixed one is favored according 
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to the MCMC fit, for botli data sets. Tlius tlie brancliing ratio to r+r~ is expected to be 
larger, wliile is suppressed. Tliis is extremely evident for Data set 1, where ~ 1/3 to e+e~ 
channel is enough to fit the data and the best-fit 5^ becomes > 0.5 (0.639 for annihilation 
and 0.655 for decay modes). Since in this case the best-fit background spectrum is about 
3.15 after propagation, compared with the one ~ 3.35 of the fixed background, the required 
electron contribution from DM is much softer. However, a non-negligible branching ratio to 
e~^e~ is still needed to describe the peak of ATIC data. A combination of fi~^fi~ and t~^t~ is 
also favored for Data set II, but with r+r" channel dominating instead, compared with the 
fixed background approach. For other parameters the quanlitative constraints are almost 
unchanged for different background choices. 

The changes of the fitting branching ratios between different background choices suggest 
that we should be cautious when claiming the DM properties used to explain the data, given 
the uncertain knowledge about the background. While global fitting will be the only way to 
seperate the background and possible "signal" , and derive model-independent implications 
of the data. 




2.2 2.3 2.4 2.5 0.4 0.6 0.8 1 1.2 2.4 2.45 2.5 0.9 1 1.1 1.2 

b - G - b - c + 



FIG. 5: Probability distributions of the eleven parameters in annihilation DM scenario used to fit 
Data set I {left) and II {right) respectively, for varying e+e~ background approach. 

With including additional parameters, we find that the constraints on each parameter are 
broader comparing with the fixed backgrond approach, however, the goodness of fit defined 
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FIG. 6: The same as Figj5]but for decaying DM scenario. 

TABLE III: Mean 1 a errors or 95% limits for the parameters in annihilation DM scenario, for the 
varying e^e~ background approach. 



Parameters 


Data Set I 


Data Set II 


my, (TeV) 


n 79Q+0.063 
U. (ZO_Q 


9 99-1+0.141 
^•^^i„0.184 


{av) (10-23 cm^ s"^) 


^■^^"-0.311 


15.625_-[^ 972 


Be 


n oc:c:+0.029 
U.OOO_o,034 


< 0.017 




< 0.486 


Q o.(,+0.027 


Br 


0.459^^^^ 


U.DD/_o.029 


7 


< 0.472 


< 0.402 


(10"5 GeV-i cm-^ s"! sr"!) 


9.793tJ:l| 


< 16.565 


6^ 


9 97n+0-036 

Z.Z/U_g 034 


9 9-1 1+0.057 
Z.Zii_o 033 


a- (10-26 GeV-i cm^^ g-i) 


1 Qqn+0.095 


2.609l°i?^ 


K 


9 9Q7+0.052 
Z.zy /_g 


2.445 ± 0.020 


c+ 




1.082 ±0.063 



as G.O.F = x^/d.o.f is getting better (TablVl). 
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TABLE IV: Mean 1 a errors or 95% limits for the parameters in decaying DM scenario for the 
varying e~^e~ background approach. 



Parameters 


Data Set 1 


Data Set II 


my (TeV) 


1.420+°-°^I 


4.501+°-305 

— U.oUo 


r (10^6 s) 


1.198+°-}0^ 

— U. iUb 


0.483+n42 


Be 


n 010+0.026 
— U.UoU 


< 0.026 


Bu 


< 0.366 


0.269+°-°27 

— U.UoU 


Br 


0.534+°-0^« 

— U.Uol 


0.715+°-°^i 

—U.UoU 


7 


< 0.917 


< 0.893 


(10-5 GeV-i cm-2 s"! si~^) 


10.924li:f,6 


< 16.377 




n r,o7+0.033 

Z.ZO(_Q 


9 O1R+0.057 


a- (10-26 GeV-i cm^^ g-^) 


1 oc;9+0.093 


9 790+0.100 


K 


9 -^ni +0.057 

Z.OUl_Q Q5g 


2.455 ±0.021 


4 


0.708l°;|20 


1.134 ±0.065 



TABLE V: The reduced Xr of each case. 





Fixed 


Varying 




Anni 


Decay 


Anni 


Decay 


Data Set I 


1.439 


1.434 


1.065 


1.075 


Data Set II 


1.143 


1.151 


1.057 


1.112 



V. SUMMARY 



In this work we employ a MCMC technique to determine the parameters of DM sce- 
nario in interpreting the recent observations of electron spectra and positron fraction data. 
Both the DM annihilation and decay are considered in this work. The DM particle is as- 
sumed to couple only with leptons. Considering the discrepancy between experimental data 
from ATIC and Fermi-LAT/H.E.S.S., we fit two data sets independently. We find that for 
Data set I (PAMELA±ATIC), DM with ^ 0.7 TeV for annihilation (or 1.4 TeV for 
decay) and a non-negligible e^e~ component is favored. For Data set II (PAMELA±Fermi- 
LAT±H.E.S.S.) f« 2.2 TeV for annihilation (or 4.5 TeV for decay) and the combination 
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of fi'^fi~ and T~^T~ final states can best fit the data. There is degeneracy between parameters 
and B^-. The constraint on the central density profile of DM halo from the diffuse 7-rays 
near the GC is also included in the MCMC code self-consistently. We find the H.E.S.S. 
observations of the GC 7-rays can give strong constraint on the DM density slope of the 
central cusp for the annihilation DM scenario. In that case the NFW profile with 7 = 1, 
which is regarded as the typical predication from cold DM scenario, is excluded with a high 
significance (> 3a). For the decaying DM scenario, the constraint is much weaker. 

We find that the largest uncertainties on the determination of DM branching ratios 
come from the limited understanding of the e~^e~ background. Besides using the canonical 
background expectations from CR propagation models, we further consider the approach 
to involve the background uncertainties. Three additional free parameters are adopted to 
describe the background contribution and fitted globally with the DM parameters. It is 
shown that the branching ratios depend sensitively on the background results. The limited 
knowledge about the background actually prevents us from a precise determination of the 
DM final state information. However, with the accumulation of more observational data 
from PAMELA, AMS02 and many other upcoming experiments with much higher precision 
and detailed information of each species, we may expect to get better understanding on both 
the background and the "exotic signal" (if exists). While global fitting will be the only way 
to separate the background and signal, and derive model-independent implications of the 
data. Since the MCMC method works in a full high-dimensional parameter space, the code 
developed in this work will be a useful tool in studying both the CR physics and the exotic 
new physics, given more precise observational data in the future. 

We do not include the radio constraints on the synchrotron radiation of the models in this 
work, though the radio observations in the GC region might set even stronger constraints 



than 7-rays, as studied in Refs. 43|-|45| . The radio constraints are found to be sensitive to the 



DM profile, instead of the central most configuration of the magnetic field. However, since 
the knowledge of the magnetic field relies on some assumptions such as the accretion state 
of central black hole and the equipartition of the matter kinetic energy with the magnetic 
pressure 4^, we think there are still uncertainties of the radio constraints. Additionally, one 
expects that antiproton constraints would be important for the tau channel, which has not 
yet been considered in the current work. This combined with the radio bound may provide 
a stringent constraint for the r channel, which deserves a detailed study in future. 



18 



As an end, we would like to mention that, there exists a self-consistent way to explain the 
current data by introducing a light mediator state in DM annihilation (or decay) 46j, which 
has not yet been considered in the current work. The distinct advantages of this type of 
models are that, the strongly constrained decay into hadrons is kinematically forbidden, and 
Sommerfeld enhancement with small velocity of Galactic DM allows for the large annihilation 
cross sections required by observations. This kind of model is not included in this work due 
to the fact that, we need to introduce more parameters in the fitting such as the mass of 
the mediator which can not be effectively constrained by the current data on the one hand, 
and the present work is the first step of the application of MCMC method in DM indirect 
searches on the other hand. We hope the forthcoming improvement of the method can 
partially address this issue. 
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